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We use modern jet-substructure techniques to propose LHC searches for multijet-resonance signals 
without leptons or missing energy. We focus on three-jet resonances produced by _R-parity-violating 
decays of boosted gluinos, showing that shape analyses searching for a mass peak can probe such 
gluinos up to masses of ~ 750 GeV (650 GeV) with 20 fb~^ (5 fb"^) at the LHC at 8 TeV. This com- 
plements existing search strategies, which also include counting methods that are inherently more 
prone to systematic uncertainties. Since 7?-parity-violating gluinos lighter than all squarks hadronize 
before decaying, we introduce new color-flow variables, "radial pull" and "axis contraction", which 
are sensitive to the color structure of the i?-hadron's decay. The former measures the inward pull 
of subjets in a fat jet, while the latter quantifies the inward drift of the A'^-subjettiness axes when 
changing the distance measure. We show that they can dramatically improve the discrimination of 
a boosted gluino signal versus QCD, tt, and combinatoric background for rrig ~ mt. Cuts on axis 
contraction also noticeably improve the resonance shape for heavy gluinos with mg > 500 GeV. 
With minor adaptations, these variables could find application in substructure searches for particles 
in different color representations or with other decay topologies. 



I. INTRODUCTION 

The Large Hadron Collider (LHC) is setting ever-more 
stringent constraints on many Beyond-Standard Model 
(BSM) theories. The most constrained BSM theories 
are those that produce large amounts of missing trans- 
verse energy (MET) and/or leptons [iH^. However, 
there are many theories that do not feature these sig- 
natures. One of the more experimentally challenging 
signals are jets (TU] with no MET or leptons, for which 
the background from ordinary Quantum Chromodynam- 
ics (QCD) processes is prodigious. 

The particular signal that we study in this paper can 
be phrased in terms of a simplified model. In addition to 
the Standard Model particle content, we consider a col- 
ored particle that is pair-produced and decays to three 
light-flavored quarks [TTUlT] |109j . A useful benchmark 
model for this scenario is the Minimal Supersymmetric 
(SUSY) Standard Model with baryon-number violating 
i?-parity violation (RPV) [HI [TH] and a gluino as the 
lightest supersymmetric particle (LSP): the colored par- 
ticles are gluinos, which each decay to a quark and an 
off-shell squark that decays to two quarks via an RPV 
coupling (see Fig. [T]). This leads to a six-jet signal, with 
two three-jet resonances from the two decaying gluinos. 
For other examples of multi-jet searches, see pPH^B] . 

If leptons and neutrinos appear in RPV cascade de- 
cays, bounds on superpartner masses can still be '-^ 
1 TeV [57]. However, the bounds are weaker if the signal 
is entirely hadronic, which we consider here. In partic- 
ular, for a gluino LSP decaying to three jets, a model- 
independent bound on exotic color octets from LEP ex- 
cludes gluino masses below 51 GeV PSI, while searches at 
CDF and CMS have only excluded gluinos with masses in 
the range ~ 77-144 GeV [29], - 200-280 GeV [30], and 
~ 280-450 GeV [31] ■ A recent ATLAS search has closed 




FIG. 1: The signal studied in this paper; pair production of 
gluinos, with each decaying to three quarks via an intermedi- 
ate off-shell squark. 



the gap between 144—200 GeV and excludes gluinos up to 
660 GeV [32] . The weakened bounds relative to i?-parity- 
conserving SUSY make RPV an attractive possibility for 
natural SUSY models ^MM- 

The six-jet signal is challenging to see due to the large 
QCD background. A further difficulty is the large combi- 
natorial ambiguity in correctly identifying the three jets 
from each gluino. The method in [iT] [29H3T] uses a cor- 
relation between the sum of the transverse momentum 
{pT,jjj) and the invariant mass (Mjjj) of three jets to 
select a phase space region with a high number of signal 
relative to combinatorial and QCD background events. 
Even so S/B often ends up being very small, motivating 
the investigation of complementary analysis methods. 

The ATLAS study [31] exploits the fact that the jets 
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from gluino decays tend to have similar px, while QCD 
six-jet events exhibit a pt hierarchy. Counting the num- 
ber of events on the high-end tail of the sixth-hardest-jet- 
Pt distribution excludes RPV gluinos up to 660 GeV ^52] , 
but does require an extremely reliable understanding of 
the background normalization. It is thus desirable to 
cross-check this result using an orthogonal search channel 
with a nearly independent set of systematic uncertainties 
(as also suggested in 

This motivates the focus of our paper: we examine 
the decay and radiation pattern produced by two boosted 
gluinos, whose decay products tend to be coUimated and 
fall into the same region of the detector, producing two 
hard fat jets. The study of boosted gluinos allows for 
a reconstruction oi rrig in the fat -jet invariant mass dis- 
tribution, giving a cleaner and more robust signal than 
other methods, and allowing a shape analysis to extract 
the gluino mass peak. We study boosted gluinos using 
jet-substructure techniques |37H61j . which have matured 
enormously in recent years and are being verified ex- 
perimentally [U [21 m |6l [62H70] (for some recent reviews 
see [7Iti75] V We achieve the best signal sensitivity using 
7V-subjettiness [371 [SH] to isolate three-pronged fat jets, 
requiring two high-px fat jets with similar masses, and 
vetoing events with a large sub-jet px-hierarchy. In ad- 
dition, we introduce two new color-flow variables, radial 
pull and axis contraction. Radial pull is based on the pull 
variable defined in |39j . while axis contraction exploits 
the shift in the minimizing axes of A'^-subjettiness when 
changing the distance measure. Both variables help dis- 
tinguish the signal's QCD radiation pattern when com- 
pared with the QCD background. 

The ATLAS study [32] examined boosted gluinos and 
excluded nig ^ 100 — 300 GeV (light compared to v^), 
but we show that a boosted search is viable at the 8 TeV 
LHC (LHC8) up to rUg < 750 GeV, where a~gg ~ 0.2 pb 
and the boosted fraction is only ©(few %). This shows 
that looking for relatively heavy gluinos in the boosted 
regime, first proposed for the Tevatron [74j, carries over 
to the LHC in spite of the smaller boosted fraction and 
production cross section of a pp collider relative to a pp 
collider. We also define a search for top-mass gluinos 
(rrig ~ TOt) with spectacular background discrimination 
that would improve on the ATLAS limit. This is relevant 
for other simplified models. 

The new variables, radial pull and axis contraction, are 
designed to measure the distribution of the soft QCD ra- 
diation pattern inside each boosted-gluino fat jet. They 
should generalize to other examples of boosted jet stud- 
ies [75], for example hadronic RPV decays of neutralino 
LSPs (as in |76j). Their use in our study relies on the 
fact that gluinos decaying via an off-shell squark gener- 
ically live longer than the hadronization time-scale, so 
that they form a color-singlet i?-hadron before decay- 
ing. This leaves a measurable soft QCD radiation pattern 
within each boosted fat jet, which differs from a beam- 
connected color octet like that for an un-hadronized 
gluino, combinatorics background, a boosted hadronic 
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FIG. 2: (a) The color flow in i?-meson (gqq) decay due to 
RPV gluino decay via an intermediate oflF-shell squark. (b) 
and (c): The two physically distinct possibilities for color flow 
in ii-baryon (gqqq) decay. Forward/backward arrows indicate 
color/anticolor, which are also indicated with r,g, b. 



top quark, or QCD background. This effect of color- 
connected jets has been previously studied for a color 
singlet (such as W or h) decaying into two jets forming 
a color dipole [SHHUl [77H5U] , and Tevatron experimen- 
tal results demonstrated its viability [52]; however, color 
fiow has to our knowledge never been investigated for 
decaying i?-hadrons or non-dipole configurations. 

In f|lll we discuss the color flow in the gluino i?-hadron 
decay and introduce variables sensitive to the radiation 
pattern. {Ill describes our Monte Carlo (MC) genera- 
tion. §IV| details the substructure variables that distin- 
guish boosted gluino decays from background, and show 
the results for searches for heavy {nig > 500 GeV) and 
top-mass gluinos. ^V| contains our conclusions. A brief 
review of RPV gluinos is in Appenix |A} 



II. NOVEL PROBES OF COLOR FLOW 

Gluinos that decay to three quarks through an off-shell 
squark and the baryon-violating RPV coupling (which 
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FIG. 3: Exaggerated schematic representation of the hypothesized radiation pattern in the plane produced by the _R-hadron 
decays shown in Fig. [2] The first radiation pattern from _R-hadron formation with its resolved color connection to the beam is 
shown in purple at the center. The centers of the hard jets from RPV giuino decay are shown in red, with orange indicating 
the shape of the soft radiation pattern resulting from the hadronization of the three quarks produced in the giuino decay. 
Squares and triangles indicate the location of axes minimizing the N-subjettiness variables t^^^ and t!^^'^ , respectively, and 
the cross marks the fat-jet centroid, i.e. the original i?-hadron direction. There are two possibilities for the radiation pattern: 
the left plot corresponds to the i?-meson and the i?-baryon in Fig.[2|a) and (b), respectively; the right plot only occurs for the 
i?-baryon in Fig. [2|c) and is of limited importance, since i?-baryons make up only ~ 1% of produced _R-hadrons [81j . 



appears in the superpotential as T4^ D X"udd) generi- 
cally have a Ufetime that is prompt but longer than the 
hadronization scale — see review in Appendix |A (We do 
not consider giuino decays originating from a displaced 
vertex or occurring outside the detector.) This means 
they first form color-singlet i?-hadrons before decaying 
to three jets. The pattern produced by the radiation 
from these signal jets is different from QCD background 
jets or from jets that originate from a particle that is 
a color fundamental (like a top quark) or octet (like an 
un-hadronized giuino decaying to three jets). In this sec- 
tion, we first give an intuitive explanation for the radia- 
tion pattern before introducing two new variables, radial 
pull and axis contraction, that attempt to quantify this. 
The new variables should, with slight adaptations, also 
be useful to distinguish background jets from jets orig- 
inating from other particles like the Higgs or W^-boson, 
which we study elsewhere |75j . 



A. Color Connections in 7?-hadron Decays 

The soft radiation pattern of a jet produced during 
hadronization will depend on how the color of the par- 
ent quark or gluon is connected to the color of the 
other quarks or gluons in the event. Put simply, the 
radiation pattern of a jet will be, on average, slightly 
'pulled' towards other jets (or the beam) to which it is 
color-connected. Furthermore, gluons have a larger color 
charge than quarks and will produce on average a some- 
what wider radiation pattern. A number of variables 
make use of this to, for example, tag dipole-pairs of jets 



and distinguish quarks from gluons [3i 

The color flow of giuino pair production is identical to 
that of gluon pair production, with gluinos being color- 
connected to the proton remnants in the beam. There- 
fore, if the gluinos were stable, we would expect the R- 
hadrons to be surrounded by a soft radiation pattern 
not unlike that of gluon jets. However, for decaying R- 
hadrons the situation is more interesting. Fig. [2] shows 
the color flow of a decaying i?-meson gqq and i?-baryon 
gqqq. Of the three hard jets from the decaying i?-meson, 
one forms a 'mesonic' singlet with one of the i?-meson 
remnant spectator quarks, while two form a 'baryonic' 
singlet with the other one of remnant quarks (a similar 
pattern of color connection with the remnants holds for 
decaying gluinoballs gg). For i?-baryons, there are two 
distinct possibilities: two baryonic singlets, one formed 
out of one hard jet and two spectators and the other 
formed out of two hard jets and one spectator; or three 
mesonic singlets incorporating one hard quark each. 

The resulting soft radiation pattern expected in the 
decay of a boosted i?-hadron is shown in strongly exag- 
gerated schematic form in Fig. [Sj The soft radiation field 
(purple region) of the first round of hadronization (dur- 
ing _R-hadron formation) is contained in the center of the 
fat jet, and might contain evidence of the initial gluino's 
color connection to the beam, indicated here by its left- 
ward slant. However, this radiation field is extremely 
soft, as well as subject to the possible shortcomings of 
i?-hadronization models, and it will not be our focus. 
Superimposed on this initial radiation field is the result 
of the i?-hadron decay, and a second round of hadroniza- 
tion to resolve the color connections amongst its decay 
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FIG. 4: Left: Normalized distributions of the new color-flow variable radial pull {t^, Eq. ([2|) for a gluino _R-liadron with 
nig — 175 GeV decaying to three jets via RPV (red and orange); QCD (green) and ti (cyan) backgrounds are also shown. The 
signal is divided into two components: red (orange) corresponds to events in which the hardest two fat jets are aligned (are 
not aligned) within AR = 0.3 of the two gluino ii-hadrons at the MC truth level. The red distribution is dominated by events 
where both fat jets reconstruct a decaying gluino, and thus constitute "good" (T) signal events, whereas the orange distribution 
shows events where the fat jets do not reconstruct a decaying gluino and are "bad" (F) signal events. For each event with 
two hard fat jets, the larger of the two radial pulls is shown in the histogram. The inset table shows the absolute sizes of the 
different samples, normalized to the number of expected events at LHC8 with 20 fb~^. Basic generator- and trigger-level cuts 



have been included, same as for Fig. 10 (a). Error bars indicate MC statistical uncertainty. 



Right: Same as plot on left, but now showing the normalized distribution of the new color-flow variable axis-contraction {A^ 
with P = 2, P' = 1, N = 3) for mg — 650 GeV. An equivalent set of cuts has been applied, same as for Fig. [7] (a). In both cases 
the distribution of the good signal events differs markedly from the other (background) distributions. 



products. It is in the radiation fields of the three hard 
subjets of the fat jet that it may be possible to find evi- 
dence of the overall singlet nature of the i?-hadron: there 
should be an overall pull towards the center, and possi- 
bly between two fat jets that form a baryonic singlet with 
part of the remnant. While there is substantial overlap 
in the signal and background distributions of color-flow 
variables, we do find that aggressive cuts on color-flow 
variables that isolate events with the patterns shown in 
Fig. [3] are very helpful in improving signal-to-background 
ratios for RPV gluino searches. Any such variable can 
also be viewed as a generalization of existing dipole tag- 
gers, and could also find application elsewhere. 



B. Radial Pull 

The distribution of the radiation field of a jet can be 
measured with a variable called pull |39j , defined as 

^ Ft 

where r is a vector in the rj — (f) plane pointing from 
the jet axis to the i*^ jet constituent, and p!^^ {p\) is 
the transverse momentum of the jet (i**^ jet constituent). 
The direction of the pull vector indicates an overall slant 
in the jet's radiation distribution, while the magnitude 
contains limited information |39| . 

For a fat jet with N subjets, we can calculate the pull 
vector for each subjet and combine them in a quantity 



that we call radial pull, defined as 
1 ^ - 

= — ij ■ fij , (radial pull) (2) 

where tj is the pull vector (normalized to unity) for the j- 
th subjet and fij is the unit vector from the fat-jet center 
to the subjet axis, i.e. the difference between the subjet's 
and fat-jet three-momentum vector, normalized to unity. 
By definition, <r G [^1) !]• 

Radial pull characterizes the degrees to which a fat jet's 
subjets are color-connected to the fat-jet center. Radial 
pull is expected to be closer to — 1 for the decay of a color- 
singlet i?-hadron located at the fat-jet center (see Fig.|3|, 
whereas fat jets from QCD are expected to have compo- 
nents with color connections to the beam and other jets 
in the event, yielding a radial pull closer to -1-1. This 
is confirmed by MC simulations of the signal and back- 
ground, which are described in detail in j |III[ Fig. |4] (left) 
shows the normalized distribution of for a gluino with 
rUg = 175 GeV (red and orange) and QCD (green) and 
ti (cyan) backgrounds (basic generator and trigger level 
cuts have been included). We divide the signal into two 
components: red (orange) in Fig. |4] corresponds to those 
events in which the hardest two fat jets line up (do not 
line up) within Ai? — 0.3 with the three-momentum of 
the two gluino i?-hadrons at the MC truth level. The 
red distribution is thus enriched with events where two 
fat jets originate from a decaying gluino, and thus con- 
stitute "good" (T) signal events. We see that the good 
signal events peak at t,- values closer to —1 and have a 
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distinct distribution from the "bad" (F) signal events and 
from QCD and ti backgrounds. 

We also tried an alternative definition of radial pull, 



N 



N 



(radial pull, alternative def.) (3) 



This variable is most sensitive to the dominant sub jet 
pull, and so i,. is peaked strongly at ±1. We find that 
it works less well in distinguishing signal from back- 
ground compared to the definition in Eq. ([2]). How- 
ever, Eq. ([3]) is better at distinguishing radiation patterns 
from hadronized gluinos (i?-hadrons) vs. unhadronized 
gluinos, see Fig. [5] It is possible that for other types of 
sig nal Eq. ^ may be more useful, but in this paper, we 
use Eq. 



C. Axis Contraction 

A second variable to measure the QCD radiation pat- 
tern of a boosted heavy particle exploits previously un- 
explored properties of A'^-subjcttiness, [37 , which 
we first review, is the average of the px's of jet- 
constituent particles, weighted by the distance to a set 
of N axes, 



^ J2 PT.min [(Ai?i,,)^ . . . , {ARn,0^] , (4) 

(5) 



■i 

do = ^ PTiRq, 



where Rq is the fat-jet radius and the sum is over all jet 
constituents i. characterizes how well the radiation 
in the jet is aligned along these TV chosen axes, which are 
labelled by a — 1, . . . , N . ARa.i is the distance in the 
rj — (j) plane between axis a and constituent i, and the 
axes are chosen in each instance to minimize the value of 
Tn I^^I unless indicated otherwise. (This is implemented 
in the Fast Jet iV-subjettiness plug-in and provides 
substantially better discriminating power in boosted par- 
ticle taggers than when the iV-subjettiness axes are cho- 
sen using a traditional jet clustering algorithm such as 
anti-fcx) j37j.) /3 determines the weighting in the sum of 
radiation far from the axes; typically /? ^ 1 — 2. 

By comparing values of r^v for different N, it is pos- 
sible to characterize the number of subjets in the event. 
For example, if T21 = tiIt\ <C 1, the radiation in the 
jet is clustered around two separate axes within the 
jet, implying that the jet contains two prominent sub- 
jets. By contrast, if T21 ^ 1, then the radiation is dis- 
tributed fairly evenly around the central jet axis and is 
not well characterized by two subjets. The ratio T21 is 
useful for tagging jets from the decay of a boosted W- 
boson, whereas the ratio T32 = t^/t2 can isolate jets with 
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FIG. 5: Comparing the normalized distributions of max(ti) 
(maximized over the two hardest fat jets in the event) for 175 
GeV gluinos with and without formation of i?-hadrons before 
decay. The kinematic cuts (though not the radial-pull cut) 
of the top-mass gluino analysis outlined in Section [IV D| have 
been applied. Comparing the first and last bin occupations 
readily distinguishes between gluinos that either form or do 
not form i?-hadrons before decaying. 



three distinct subjets, which is useful to tag boosted top 
quarks [37l |38l [HI [7l] as well as boosted RPV gluino 
decays, as we show below. 

Varying the parameter /3 changes the sensitivity of 
to radiation far from the TV-sub jettiness axes, and can 
therefore probe the shape of the radiation inside the fat 
jet. This provides a sensitivity not just to the num- 
ber of 'hard subjet-like structures' in a fat jet (which 
might be accessed by reclustering) but also to the shape 
of their radiation patterns. Traditionally, this has been 
exploited in top-taggers [37j[38] by setting /3 w 1 to make 
iV-subjettiness hone in on very tight radiation centers. 

A novel way of exploiting this shape sensitivity is to 
study the j3 -variance of the axes which minimize t^. 
(The idea of varying the parameters of event shape vari- 
ables to extract additional information was first proposed 
in [55].) Consider a fat jet with N well-defined subjets. 
We denote the a = 1, . . . , TV axes which minimize by 



R 



= (<7V>f,Ar)- 



(6) 



If the soft radiation field of a subjet associated with axis 
a is skewed towards one side, then the axis is shifted away 
from the center of the subjet, and this shift is more pro- 
nounced for higher values of /?. Therefore, for example, 
the vector 



AR' 



■a, AT 



R 



a,N 



R 



0' 

'a,N' 



{axis pull) (7) 



which we call axis pull, should point in the direction of 
the skew for /3 > /3' . This typically points in the same 
direction as the pull ([I]). 

For the radiation pattern from gluinos, shown in Fig.jSj 
the skew of the radiation field should lie between the 
subjet axis and the fat jet center. To capture this effect. 
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we define a scalar quantity called axis contraction as 



— ^ 



N 

E 



R 



a.N 



= 1 

E| 

a=l 



R 



0' 

a.N 



Rc 



Rr 



(axis contraction) 



(8) 

where i?ccn rnarks the fat jet centroid, i.e. the fat jet 
momentum in the rj — cj) plane, and j3 > P' . For our 
application we choose (/3, /3') = (2, 1). 

If the axes shift towards (away) from the centroid of 

the fat jet as /? is changed from 1 to 2, then A^^ is smaller 
(larger) than 0. We therefore expect that for boosted 
RPV gluinos, is on average smaller for signal fat jets 
that contain all of a gluino's decay products than for fat 
jets from kinematically identical QCD or tt backgrounds, 
or signal fat jets that do not contain the decay products 
of a single gluino. This is again confirmed by MC simula- 
tions of the signal and background, which are described in 
detail in { III Fig. |4] (right) shows the normalized distri- 
butions of A-^ for a gluino with mg = 650 GeV (red and 
orange) and QCD (green) and tt (cyan) backgrounds (ba- 
sic generator and trigger level cuts have been included). 

Note that the small numerical range of the axis con- 
traction variable over our event samples is not indicative 
of any unrealistic detector resolution required to observe 
its variation, but is just a matter of how the variable is 
normalized. Furthermore, the shapes of the jet-pT dis- 
tributions are very different at maximal vs minimal val- 
ues of axis contraction (or radial pull) . These differences 
are easily distinguished by the calorimeter and particle 
tracker. 

The idea of using a shift in the A^-subjettiness axes un- 
der a change of j3 to probe color flow is very general. One 
could imagine using A^ as a dipole tagger, for example. 
This will be investigated further in |75| . 



III. MONTE CARLO EVENT GENERATION 

We simulated gluino pair production, i?-hadron for- 
mation, and RPV decay with a developmental version 
of Pythia 8.165 pti55] [TTU] . This was done for two 
ranges: "top-mass" gluinos with rUg ~ 175 GeV and 
heavy gluinos with m-g = 500, 550, . . . 1000 GeV. The 
number of unweighted signal events simulated for the 
heavy (light) gluino analysis was ~ 10^ — 10^ (^ 4 x 10^) 
per mass point. All distributions were reweighted to cor- 
respond to a given luminosity at LHC8, usually 20 fb^^. 
Samples without i?-hadron formation before the decay 
were also generated for comparison. All signal cross sec- 
tions were calculated at NLO in Prospino [87] . 

We used Sherpa 1.4.0 [HI |89j to generate fully 
matched QCD (with 2-6 hard jets from the Leading- 
Order matrix element) and ti background samples (fully 
hadronic decays with < 2 additional hard jets), with ad- 



ditional jets produced by the shower. Event generation 
was weighted to adequately sample high-multiplicity and 
high-pT events. We generated separate backgrounds for 
the heavy and the top-mass gluino searches: 

• Heavy gluino (rUg > 500 GeV) search: The relevant 
trigger is Ht ^ 850 GeV along with at least one 
jet of pt > 200 GeV [90], while the analysis itself 
requires two hard fat jets. Accordingly, we required 
background events at generator level to have two 
fat jets (anti-fcx, R = 1-5) each with pT > 500 
GeV, along with one thin jet (anti-fcx, R = 0.4) 
with pt > 200 GeV. All generator-level cuts act 
at parton-level, and thresholds are correspondingly 
conservative compared to later cuts on the fully 
showered event. Because efficiencies of subsequent 
cuts are low, we generated a sample of ~ 50 million 
QCD events and ~ 10 million tt events. We find 
that the dominant background by far is QCD. Top 
backgrounds have a much smaller cross section and 
tend to produce fat jets that lie around ^ 175 
GeV, well below the gluino mass scale under study. 

• Top-mass gluino search: This search uses a trig- 
ger requiring six thin jets (anti-fcT, R = 0.4) with 
Pt ^ 60 GeV for full efficiency. Therefore back- 
ground events were required at generator level to 
have two fat jets (anti-Zcxi R — 1-5) with > 150 
GeV, as well as at least four thin jets with p^ > 40 
GeV. Nevertheless, the poor trigger and cut effi- 
ciencies required the generation of around 100 mil- 
lion QCD events. The top, for which we generated 
~ 20 million events, is subdominant to QCD for 
most of our subsequent chain of cuts, but can dom- 
inate after the final color-flow cuts. 

All events are analyzed in our own Fast Jet 3.0.2 
91|-based code, interfaced with the A^-subjettiness plu- 
gin [38] to calculate and find minimizing axes. His- 
tograms were combined and reweighted in Mathematica 
to the number of expected events at LIIC8 for a given 
integrated luminosity. For a histogram bin containing 
events with weight {wi}, the MC statistical error on the 
sum Wi is wf, which we include in our distribu- 

tions. We generated enough events to ensure that this 
statistical error is no larger than Poisson error corre- 
sponding to the actual expected number of LHC8 events, 
and that the distributions are as smooth as possible. Sys- 
tematic errors are not included in our predictions, since 
the number of events surviving all the cuts is generally 
so small that Poisson uncertainties dominate. Detector 
effects are not explicitly modeled, but are unlikely to be 
a limiting factor since we only use the central part of 
the event {\ri\ < 2.5), our fat jets are extremely hard 
{pt > 600 GeV) and low statistics necessitate a large 
bin size for fat-jet-mass histograms anyway, precluding 
the use of any unrealistic mass resolution. Furthermore, 
as mentioned in the previous section, our color-flow vari- 
ables do not rely on any unrealistic detector resolution 
to achieve their distinguishing power. 
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itly demonstrate that our cuts eliminate the QCD and tt 
as well as the combinatorics background. 



IV. ANALYSIS 

This section details our analysis strategy and results. 



A. Strategy 

We discuss here the set of variables used to distinguish 
signal from background. Our aim is to remove as much 
QCD background as possible and have the gluino reso- 
nance be clearly distinguishable from the background. In 
particular, we do not rely on a counting analysis, as this 
is more susceptible to uncertainties in the background 
normalization. 

In addition to the generator- and trigger-level cuts dis- 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 cussed in § |III[ we select events with the following prop- 

T32 erties: 



FIG. 6; Comparison of 300 GeV< pt < 400 GeV, R = l.O 
anti-fcT fat jets from Sherpa and the ATLAS 7 TeV, 35 pb~^ 
analysis [63j- The plots are normalized to the same area. 
Top: Jet-mass distribution. Bottom: A''-subjettiness ratio 
T32 = T3/T2 distribution. 



All signal and background samples include the under- 
lying event and are generated using default parton dis- 
tribution functions. To verify their validity, we applied 
several cross-checks to our background samples. We were 
able to reproduce the QCD 6-jet distributions from [SO] 
as well as jet-mass and iV-subjettiness distributions from 
[55] (see Fig. |6]). While the shape matched well across 
the entire range of all distributions (particularly for A^- 
subjettiness), the normalization was half as large as the 
data. Therefore, we scaled up the cross section of our en- 
tire QCD background sample (supplied at LO by Sherpa) 
by a factor of 2. It was unnecessary to similarly re- 
scale the top background: it is negligible for the heavy 
gluino search, and while it is the dominant background 
after color cuts are applied in the top-mass gluino search, 
it is nevertheless eliminated by an additional 5-veto. 

Only ^ 10% of top-mass gluinos and ^ 1 — 4% of heavy 
gluinos with a mass of 500—1000 GeV are boosted enough 
for all the decay products to potentially end up in two 
fat jets. Since our search attempts to reconstruct a res- 
onance in the fat jet mass spectrum, the large major- 
ity of signal events has to be treated as combinatorics 
background. Therefore, we divide the signal sample into 
two groups: "good" (T) signal, which is enriched with 
boosted gluinos by requiring that the two hardest fat 
jets are within AR = 0.3 of the truth-level gluinos, and 
"bad" (F) signal, i.e. the rest. This allows us to explic- 



1. Two hard fat jets of size AR — 1.5, clustered with 
the anti-fcr algorithm (92] • Each fat jet contains 
the three jets from a gluino decay if the gluino is 
boosted enough. We always consider only the two 
hardest fat jets in each event, and in many cases cut 
on the larger/smaller value of a variable evaluated 
for each jet (indicated with max/min). 

2. Each fat jet should be three-pronged. We use the 



N-subjettiness variable defined in Eq. (|4| 
lect events with small values of 



to se- 



T32 



(9) 



3. The two signal fat jets should have comparable in- 
variant masses, since both reconstruct the gluino 
resonance. We thus select events with small values 
of the jet-mass symmetry parameter 



\mi - m2\ 
{nil + m2)/2 



(10) 



4. Each fat jet should contain three hard subjets with 
similar px, since the QCD background jets tend to 
have asymmetric pt^s. We thus select events where 
the subjet hierarchy parameter 



^31 = 



PT3 
PTl 



(11) 



is close to one. Here pTi is the i-th subjet found 
by reclustering each fat jet constituents with anti- 
kr, R — 0.4. This variable is similar to the flow 
variable in [93) . Versions of both and /131 have 
been used in 



8 



5. The radiation pattern within each fat jet should be 
that of a color-singlet _R-hadron. We thus select 
events with constraints on the new variables radial 
pull or axis contraction defined in Eqs. ([2| and (|8]), 
respectively. For radial pull, we use subjcts ob- 
tained by reclustering the fat-jet constituents with 
anti-kT, R = 0.4, the same as for h^i. 

We also evaluated other substructure variables like 
girth [IS |44] and planar flow [45l [94] , but they did not 
add any discriminating power in our case. 



B. Pile- up Considerations & Effect of 
Jet-Grooming 

In the high luminosity environment of the LHC, the 
effects of pile-up (PU) can distort the spectrum of jet 
observables such as its mass and iV-subjettiness. The 
effects of pile-up can be largely eliminated using one of 
several jet-grooming techniques [1S1 - I5U1 . While we do 
not explicitly include the effects of jet-grooming in our 
main analysis, we verify here that an analysis based on 
variables outlined in pV A| givc results that are applicable 
when jet grooming is included. 

To see the effects of jet grooming, we trim [33] our 
samples by re-clustering fat jets into subjets of i? = 0.4 
using the kj^ algorithm and discarding all subjets with 
PTsubjet/pTjct < /cut- We take /cut = 0.02 (/cut = 0.05) 
for the heavy (top-mass) gluino analysis, with the precise 
value chosen to eliminate most of the effects of PU based 
on a simple model with mean PU energy of 12 GeV per 
unit area and intra-event fluctuations of 3 GeV per unit 
area (characteristic of PU with 20 primary vertices per 
event) [96]. We choose a smaller /cut for heavier gluinos 
than for rrig ~ rut, since PU remains fixed even as the 
fat-jet pt increases for heavier gluinos. These parameters 
are comparable to those used by [52] , 

Trimming does alter the distributions of some of our 
kinematic variables; in particular, it tends to shift A^- 
subjettiness to lower values. However, we find that 
it is always possible to choose values of the cuts on 
T'32, Sm, and /i3i such that the overall signal and back- 
ground efficiencies are essentially the same for trimmed 
and untrimmed samples, and so the outcomes of the un- 
groomed analyses remain the same when grooming is in- 
cluded. In fact, it is possible in some instances to achieve 
a background acceptance that is lower for the trimmed 
sample while maintaining the same signal efficiency; this 
makes our ungroomed analysis somewhat conservative, 
although we leave an optimization of jet grooming pa- 
rameters to future work. 

Our color-flow variables, however, must be calculated 
on the jets prior to grooming, because grooming can re- 
move the soft radiation within and between subjets that 
distinguishes different color flows. To mitigate the effects 
of PU, it is preferable to include soft radiation over the 
smallest possible area. We find that, in most events, the 



values of the color-flow variables are dominated by ra- 
diation close to the subjet center, and we suggest that 
including only soft radiation within R = 0.4 of each sub- 
jet axis would reduce contamination from PU (this is 
similar to the method of PU mitigation used in [^). 
Furthermore, we find that the scale of radiation driving 
the color-flow variables is ^ 12 — 30 GeV (and sometimes 
higher for heavy gluinos), which is above the character- 
istic scale of intra-event fluctuations in PU (~ 3 GeV for 
20 primary vertices [96]). We do not, however, include 
PU in our MC samples, and we leave it to the experi- 
mental collaborations to study precisely how color-flow 
variables are affected by PU and to conflrm their utility. 



C. Results for nig > 500 GeV 

Table |l] details our cuts, together with the signal and 
background efficiencies for nig G [500, 850] GeV and the 
expected number of events at LHC8 at 20 fb~^. We note, 
firstly, that due to the generator-level cuts outlined in 
Section |III[ the event numbers for background samples 
can only be compared to actual data numbers after the 
third cut, which requires two fat jets with pT > 600 GeV. 
Secondly, even though we classify the entire gluino sam- 
ple as 'signal', only the C(few %) boosted fraction that 
ends up in a fat jet that reconstructs the gluino mass 
is of interest. Therefore, while the signal efficiency is ex- 
tremely low, a more relevant figure-of-merit is the ~ 10% 
final acceptance of that boosted fraction. The cuts are 
extremely effective at reducing background from ~ 10^ 
events after the first fat-jet pr cut to ~ 10 events after 
all other cuts, and we can reconstruct gluino resonances 
as heavy as ^ 750 GeV. Finally, the last two cuts have 
been optimized for each Wg, mostly to target the boosted 
fraction (px ^ rUg) and make the cut on axis contraction 
Eq. (|8|) more conservative for larger where there is 
less (boosted) signal available. Interestingly, for the con- 
servative cuts necessary to retain high signal efficiency 
in the heavy gluino case, axis contraction turns out to 
be slightly better-suited than radial pull, which is why 
we use the former in spite of the latter's superior signal 
separation on the tail of the distribution. 

Fig. [7] illustrates the cut chain following the first fat- 
jet pt and T32 cuts, for rUg — 650 GeV. As can be 
seen from the jet-mass distribution (a), the QCD back- 
ground (green) completely dominates over both top back- 
ground (cyan) and signal (split into 'good' signal in red, 
where the hardest two fat jets are within AR = 0.3 of 
the gluinos, and the rest in orange). Each unstacked 
histogram shows the unity-normalized distributions of a 
variable before cutting on it, for each signal and back- 
ground component. The cut is indicated with a dashed 
magenta line, with events in the direction of the arrow 
kept for the next cut. (f) and (h) show the fat-jet-mass 
distribution before and after a conservative cut on axis 
contraction, shown in (g), is applied. (For a distribution 
of axis contraction at the same stage of cuts as (a), (b). 
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1.1 (0.62) 
18 (49) % 
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TABLE I: Number of expected events at LHC8 with 20 fb ^ for signal and background after each cut for the heavy gluino 



analyses, and gluino pair production cross sections [83 in square brackets. The generator-level cuts outlined in Section III 



only applied to the background samples, "max" and "min" in the cut variables applies to the two values obtained for the two 
hardest fat jets in each event. The second line in each row shows the efficiency of each individual cut step. For signal, the 
numbers in brackets refer to the "good" combinatorially correct signal defined in Section [Till The last pr and axis-contraction 
cuts were optimized for each mg, with the following thresholds {p??"^, ^'""^j: rrig = 500: {600,0.02}, rrig = 550: {600,0.02}, 
m-g = 600: {600,0.02}, m~g = 650: {700,0.04}, nig = 700: {750,0.03}, mg = 750: {800,0.03}, rUg = 800: {850,0.04}, rrig = 850: 
{900,0.03}. (Masses and momenta in GeV.) 
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FIG. 7: Distributions for kinematic and substructure variables at various stages in our chain of cuts for mg — 650 GeV. 
Unstacked distributions are separately normalized to unity, while stacked distributions show actual number of events expected 
at LHC8 for 20 fb"\ QCD (tt) background is in green (cyan), while signal is in red (orange) for events whose fat-jet momenta 
are aligned within AR = 0.3 of the gluino f?-hadron momenta checked with MC truth. The shown distributions are (a) the 
invariant mass Mj of both hardest fat jets in each event, for events containing two fat jets with pr > 600 GeV and T32 < 0.7, 
(b) the maximum N-subjettiness T32 (Eq.((9|) of the two fat jets, with the same cuts applied, (c) the jet-mass symmetry Sm 
(Eq.(lOl) after the cut max(r32) < 0.5, (djthe minimum subjet hierarchy /131 (Eq.(lll) of the two fat jets including the cut 
Sm < 0.1, (e) the minimum fat-jet momentum p^' of the two fat jets including the cut min(/i3i) > 0.2, (f) the resulting 
Mj, (g) the color variable axis contraction max(j4|^) (Eq.|8|), and (h) the final Mj including the cut max(A|^) < 0.03. MC 
statistical error is shown in faint-colored bands with the color corresponding to the respective signal or background component. 
In Mj distributions, the gray error bands indicate the MC statistical error of the various components added in quadrature. A 
vertical magenta line indicates the applied cut, with the arrow pointing to the events that are kept. The inset shows the actual 
number of expected events at LHC8 with 20fb~^ for each signal and background component. 
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FIG. 8: From top to bottom, stacked fat-jet-mass distributions of signal and background for mg — 600, 700, 800 GeV, similar 
to Fig. [7](f) and (h). Plots on the left include the final axis contraction cut, while plots on the right do not. 



see Fig. |4] (b)). The color cut loses very little signal but 
reduces background by almost one-half, and significantly 
improves the shape of the distribution, making the gluino 
peak stand out very clearly. This is also illustrated in 
Fig. [8] for a few other gluino masses. 

To estimate our analysis' mass reach, we follow the 
maximum likelihood procedure for a shape analysis out- 
lined in [5S], but using the CLg method [HZ]. This in- 
volves using our predictions for the final jet-mass distri- 
butions, e.g. Fig. |8) to produce large collections of pseu- 
dodata for the 'background-only' and the 'background 
-|- signal' hypothesis, where the signal is scaled by some 



overall signal strength fi. (We ignore the statistical un- 
certainties of the MC prediction, as for such low event 
rates the Poisson fluctuations dominate.) By comparing 
the distributions of the CLg test statistic computed for 
those pseudodata sets one can arrive at an expected sig- 
nal strength exclusion, as well as the exclusion's ±la and 
±2(7 fluctuations. The background normalization was al- 
lowed to float in each of these fits individually, making 
this a true shape analysis. 

The resulting expected cross-section exclusions are 
compared to the gluino production cross section at LHC8 
with 20 fb~^ in Fig. [oj both with and without the axis- 
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FIG. 9: Expected 95% confidence-level cross section exclusion (and its ±1ct and ±2(t fluctuations) obtained with the CL^ 
method [97] for different gluino masses at LHC8 with 20 fb~^ of data. The left (right) exclusions were derived without (with) 
color-flow cuts on axis contraction (Eq. ([8|), with a floating background normalization in all fits to pseudodata. In both cases 
the reach is ~ 750 GeV, though the color-flow cuts improve exclusion at lower masses and improve the resonance shape. Fixing 
the background normalization increases the reach by ~ 20 GeV (not shown). For 5fb~^, the reach is ~ 650 GeV (not shown). 



Analysis 


generator 

level 

cuts 


Trigger: 

6 thin jets with 
Pt > 60 GeV 


two fat jets 
with 

PT > 200 GeV 


d 
V 

B 


d 

V 

X 

B 


d 

V 

CO 

s 


d 

A 

CO 

s 


radial pull cut 
max(tr) < -0.6 


fe-veto 


175 GeV 


QCD 


2.5 X 10» 


1.9 X 10'' 
0.78% 


1.7 x 10** 

88% 


4.1 x 10^ 

24% 


2.0 x 10* 
4.8% 


4.5 x 10^ 
24% 


488 
11% 


0.86 
0.18% 


0.81 ± 0.24 

94% 


top 


1.1 X 10'^ 


2.6 X 10* 
2.4% 


2.3 x 10* 

88% 


1.2 x 10* 
51% 


2.2 x lO'^ 
19% 


771 

36% 


253 
33% 


5.6 
2.2% 


0.48 ± 0.25 
8.6% 


II 

iCri 


Gluinos 
[1.9nb] 


3.7 X 10^ 


1.1 x 10'^ 
2.9% 


9.5 x 10^ 
90% 


3.3 x 10^ 
34% 


2.0 x 10-* 
(6.6 x 10^) 
6.1% 
(9.9%) 


6.1 x 10^ 
(3.2 X 10^) 
30% 
(49%) 


1.5 X 10^ 
(1.0 X 10^) 
24% 
(31%) 


128 
(118) 
8.7% 
(12%) 


120 ± 17 
(111 ± 15) 
94% 
(94%) 
S/B = 93 



TABLE II: Number of expected events at LHC8 with 20 fb^^ for signal and background after each performed cut for the 
top-mass gluino analysis, and gluino pair production cross section [87 in square brackets. The generator-level cuts outlined 
in Section |III| are only applied to the background samples, "max" and "min" in the cut variables applies to the two values 
obtained for the two hardest fat jet in each event. The percentage in each cell is the efliciency of that column's individual cut 
step. For signal, the numbers in brackets refer to the "good" combinatorially correct signal defined in Section|lIl] The cut chain 
after the second vertical double-line is illustrated in Fig. |10[ Note that the efficiency of the "bad" signal for the color-flow cut 
is only 2.3%, compared to 12% for the "good" signal. 



contraction cut. The expected reach is ~ 750 GeV and 
could be as high as 800 GeV, surprisingly high for a fully 
boosted search. With 5 fb~^ of LHC8 data the reach is 
rrig ~ 650 GeV, higher than the LHC7 CMS search [3T] 
and comparable to the ATLAS resolved search [32] , which 
was a pure counting experiment. Since the boosted anal- 
ysis has S/B ~ 1 and is predominantly limited by statis- 
tics, more data should lead to better mass exclusion using 
the boosted analysis techniques discussed here. 

Due to the small number of events surviving all the 
cuts, our estimate of the mass reach is not actually in- 
creased by the color-flow cuts, though the excluded cross 
section is increased for small nig, where more events sur- 
vive all cuts. Furthermore, the color-flow cuts do de- 
crease the background at little cost to the signal and 
may be useful to control systematic uncertainties on the 



background. 

Exclusion could be improved by studying the fat-jet- 
mass distribution in control samples where no signal 
events are expected; a similar procedure was applied to 
the thrcc-jct-mass distribution by CMS [3D1 [3T], effec- 
tively fixing the background normalization. This would 
increase our mass reach by ^ 20 GeV, the small improve- 
ment being indicative that our simple shape analysis al- 
ready does well in fixing the background normalization. 



D. Results for rrig ~ mt 

While the recent ATLAS search for RPV gluinos ex- 
cluded the mass range rrig ~ 140 — 200 GeV [32], we 
still consider gluinos in this mass window as they pro- 
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FIG. 10: Distributions for kinematic and substructure variables at various stages in our chain of cuts for rrig = 175 GeV. 
Unstacked distributions are separately normalized to unity, while stacked distributions show actual number of events expected 
at LHC8 for 20 fb'^ QCD (tt) background is in green (cyan), while signal is in red (orange) for events whose fat-jet momenta 
are aligned within AR — 0.3 of the gluino i?-hadron momenta checked with MC truth. The shown distributions are (a) the 
invariant mass Mj of both hardest fat jets in each event, for events passing the 6-jet (60 GeV) trigger and containing two fat 
jets with pt > 200 GeV and r32 < 0.7, (b) the maximum N-subjettiness r32 (Eq.([9|) of the two fat jets, with the sa me cuts 
applied, (c) the jet-mass symmetry Sm (Eq.(lOl) after also max(T32) < 0.5, (d) the minimum subjet hierarchy /lai (Eq.( |ll[ )) of 
the two fat jets including also Sm < 0.1, (e) the resulting Mj distribution, (f) the color variable radial pull max(fr) (Eq.([2|), 
(g) the Mj distribution after the radial pull cut, and (h) the same distribution after a 6-veto. MC statistical error is shown 
in faint-colored bands with the color corresponding to the respective signal or background component. In Mj distributions, 
the gray error bands indicate the MC statistical error of the various components added in quadrature. A vertical magenta line 
indicates the applied cut, with the arrow pointing to the events that are kept. The inset shows the actual number of expected 
events at LHC8 with 20fb~^ for each signal and background component. 
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vide a large sample of highly boosted signal events with 
which we can demonstrate the full discriminating power 
of our color-flow variables. Furthermore, we show that 
a dramatically higher signal purity than in [32' can be 
achieved, which translates to a correspondingly higher 
cross section exclusion. 

We find that the trigger requiring six jets with pT > 
60 GeV yields the highest signal efficiency and tends to 
give boosted gluinos. Due to the abundance of signal 
we can cut hard on the substructure and color-flow vari- 
ables. The cuts are outlined in Table |ll] and illustrated 
in Fig. [To] While axis contraction was better suited for 
the soft cuts on color flow required by low signal numbers 
in the heavy gluino analysis, radial pull is more effective 
when we have enough signal to cut more aggressively. 
The result is an extremely powerful color-flow cut that 
keeps 12% of the combinatorially correct signal (of events 
surviving the previous cut step) but only 2% of com- 
binatorially incorrect signal or ti background, and only 
0.2% of QCD background. The resulting mass distri- 
bution Fig. 10 'g) is very signal-dominated and displays 
a clean resonance at 175 GeV [ 111 . Applying a 6- veto 
(assuming 70% tagging efficiency for b quarks and a 1% 
light-jet mistag rate (99^) yields Fig. 10 [h), with a res- 
onance peak of ~ 100 gluino events and only about 1 
background event in total. 

It is instructive to compare the radial-pull distribution 
after the kinematic substructure cuts Fig. 10 'f) to Fig. [3] 
(left), which shows the same distribution before many of 
the kinematic and substructure cuts are applied. The 
similarity in the distributions (accounting for increased 
statistical uncertainty of the MC predictions after more 
cuts) suggests that radial pull is relatively uncorrelated 
with other substructure variables, and provides genuinely 
new discriminating power that cannot be accessed by 
simply making another cut more aggressive. 

The counting experiment performed by ATLAS at 
LHC7 with 5 fb~^ has excluded top-mass gluinos with 
Br(^ -> qqq) < 0.25. Treating our nearly background- 
free sample in a similar way would exclude Br (5 
qqq) < 0.15 at LHC8 with the same luminosity, decreas- 
ing to Br(^ qqq) < 0.05 at 20 fb~^. Apart from 
excluding other possible RPV spectra, this can be rele- 
vant for different color representations decaying to three 
jets, which have different production cross sections. The 
shape analysis serves as a useful check on any discovered 
signal, since a resonance is clearly constructed. Further- 
more, if a three-jet resonance in this mass range were to 
be found, one could use the alternative version of radial 
pull, Eq. ([3]), to distinguish a hadronized color-octet from 
some other state, as illustrated in Fig. [5] 



boosted fraction of very heavy particles, such as RPV 
gluinos. Although this has recently been demonstrated 
in an LHC search for top-mass gluinos [32]) we propose 
such search strategies also for heavier masses at LHC8. 
We use existing variables, such as iV-subjettiness, jet- 
mass symmetry, and subjet hierarchy, as well as our new 
color-flow variables, radial pull and axis contraction, to 
isolate a high-purity signal sample. Our suggested anal- 
ysis strategy with aggressive cuts on these variables can 
be competitive with existing search strategies (TT] [29l - 
I32| . while providing nearly independent systematic un- 
certainties. With 20 fb"^ (5 fb"^) at LHC8, it should be 
possible to probe boosted gluino-like resonances around 
the top mass with ct x Br < 0.3 (0.1) nb, corresponding 
to ~ 0.05 (0.15) gluino branching fractions to three jets. 
Heavy gluinos can be excluded for nig < 750 (650) GeV. 

Additionally, we have shown that radial pull and axis 
contraction appear to be powerful variables for distin- 
guishing processes with different color flows, particularly 
for highly boosted events. We anticipate that they have 
applications in studies of other boosted colored objects, 
and may prove useful in distinguishing different models 
in the case that a boosted resonance is discovered. 



Acknowledgements 

We are grateful to T. Cohen, M. Freytsis, E. Izaguirre, 
and M. Lisanti for helpful discussions and comments on 
a draft version of this paper. We also thank M. Begel, 
E. Halkiadakis, J. Hobbs, D. Tsybychev, P. Meade, 
D. Miller, G. Salam, J. Shelton, G. Sterman, J. Thaler, 
and J. Wacker for helpful conversations. We especially 
thank P. Skands and T. Sjostrand for their help with 
Pythia, and S. Hoche and S. Schumann for their help 
with Sherpa. The work of D.C. was supported in part 
by the National Science Foundation under Grant PHY- 
0969739. RE is supported by the Department of En- 
ergy Early Career research program under Award Num- 
ber DE-SC0008061. The work of B.S. was supported in 
part by the Harvard Center for the Fundamental Laws of 
Nature, the National Science Foundation under Grant 
PHY-0855591, and the Canadian Institute of Particle 
Physics. The work of D.C. and R.E. was conducted 
in part at the Aspen Center for Physics, supported by 
the National Science Foundation under Grant No. PHY- 
1066293. Some of the numerical calculations in this paper 
were performed on the Odyssey cluster supported by the 
FAS Research Group at Harvard University. 



Appendix A: Review of RPV Gluinos 



V. CONCLUSIONS 

Jet-substructure variables are excellent tools in 
searches for boosted resonances. We have shown that 
such techniques can also be useful when applied to the 



In addition to the superpotential terms of the Mini- 
mal Supersymmetric Standard Model, Standard Model 
gauge invariance allows terms that violate baryon (B) 
and lepton (L) number. B and L violation are strongly 
constrained by limits on proton decay and neutrino 



15 



masses [in]- Often, i?-parity is imposed to remove these 
dangerous terms. However, if lepton number is con- 
served, the B-violating RPV term 



(Al) 



is allowed and only loosely constrained. Many theories 
accommodate RPV [T00Hl04] . 

Models with A" ^ can be challenging to observe, 
since the hadronic decay of supersymmetric particles suf- 
fer large QCD backgrounds. For RPV operators with one 
or more heavy-flavor quarks, multiple 6-tags, leptons, and 
missing energy can help distinguish the signal from back- 
grounds [33, JSJ |35] . If the quarks are all light-flavored, 
however, no easy distinguishing property exists. 

The couplings Xljk ^re subject to various constraints 
for different {ijk}. The dominant constraint on \ii2 
comes from nucleon-antinucleon oscillation through an 
intermediate strange squark and gluino. The precise 
bound depends on the off-diagonal entries of the unknown 
left-right strange squark mixing matrix. If there is no 
mixing suppression, the bound is \\ii2\ ^ 10~® for rrig ^ 
TeV and mg^ ~ 5 TeV [15] . The bound is greatly relaxed 
if the mixing is suppressed. 

The coupling A'/ji induces strangeness-violating 
nucleon-antinucleon oscillation, leading to double nu- 
cleon decay processes such as pp — >■ K^K^ and nn — > 
K^K'^. Bounds extracted from this process are highly 
sensitive to hadronic and nuclear matrix elements, and 



range from |A'/2i| < lO"'^ to 1 [102], with ^ quoting a 
value of |A'/2i| < 10"^ for ms^ = rUg = 300 GeV. 

The couplings A212 and A221 are even less constrained. 
Bounds can be imposed by requiring that the RPV cou- 
plings remain perturbative up to the GUT scale, giving 
|A^'i2UA^'2i| < 1.25. 

We focus on models with light gluinos, which are pair- 
produced and each decay through an off-shell squark 
to three jets. Models exist in which the gauginos 
are naturally lighter than the scalars, including "split 
SUSY" [inSl UnS]), leading to the possibility of a gluino 
LSP. The collider phenomenology depends on the gluino 
lifetime. Gluinos decaying through an off-shell squark 
and RPV couplings always hadronize before decaying, and 
depending on the magnitude of A", may also give rise to 
displaced vertices. To see this, we estimate the gluino 
width in the limit rriq ^ rrig: 



|A"| 



"12 m| 



(A2) 



then Fg 



Using a conservative value |A"| = 10 ^ 
10^" GeV with m-g = I TeV and = 
Fg ~ 10~^ GeV with the most relaxed bound |A" 



5 TeV, while 
1. 

Since these are both well below Aqcdi the gluino forms 
an i?-hadron prior to decay whenever > rrig, even 
when the RPV couplings are 0(1). We consider only 
prompt decays after hadronization; for long-lived and 
stopped gluinos, see e.g. |1071I108) . 
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